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Abstract 

Using Jacobi elliptic function addition formulas and summation identities we obtain several static 
and moving periodic soliton solutions of a classical anisotropic, discrete Heisenberg spin chain with 
and without an external magnetic field. We predict the dispersion relations of these nonlinear 
excitations and contrast them with that of magnons and relate these findings to the materials 
realized by a discrete spin chain. As limiting cases, we discuss different forms of domain wall 
structures and their properties. 
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I. INTRODUCTION 



Establishing integrability and obtaining exact solutions of discrete nonlinear physical 
systems are important issues of current interest. Starting with the integrable discrete model 
of Ablowitz and Ladik for several other discrete nonlinear evolution equations exact 

n 

elliptic function and soliton solutions have been obtained in recent years [2j]. These include 
certain discrete versions of nonlinear Schrodinger (NLS) equation j^, 4|, (j)^ equation js,!^, 
derivative NLS equation [zt], coupled 0"^ equation [8i], coupled asymmetric double well and 
coupled (f)^ equation j^, complex modified Korteweg-de Vries equation [l^, etc., where 
effective use of summation relations of Jacobian elliptic functions was made and periodic 
and solitary wave solutions of moving and static types obtained. 

In this connection, a physically important discrete nonlinear dynamical system which has 
been of considerable interest in diverse areas of physics for a long time is the anisotropic 
Heisenberg ferromagnetic (and antiferromagnetic) spin system with or without an external 
magnetic field. It has been studied for various aspects in magnetism, condensed matter 
physics/materials science, statistical physics, nonlinear dynamics, etc. both from classi- 
cal and quantum points of view For example, the one-dimensional quantum spin-1/2 
XYZ chain has been shown to be an exactly solvable system either through Bethe ansatz 



procedure or through quantum inverse scattering method [12|, |13| and the eigenvalue spec- 
trum and eigenfunctions have been obtained. For large value of spins, however, a classi- 
cal/quasiclassical description has been known to be an adequate description so that spins 
can be treated as unit vectors and classical equations of motion for the spin vectors can be 
obtained as limiting forms of the quantum equation of motion or as dynamical equations 
derived from postulated spin Poisson bracket relations [l^. Another area of considerable 
physical interest in which such classical anisotropic spin systems have been studied in the 
presence of Gilbert damping is the microscopic behavior of spin waves in magnetic bodies 
of arbitrary shape 15| and the study of spin-torque effect in ferromagnetic layers with spin 
currents [Ifl] on spin waves and domain walls. Recently it has also been pointed out that 



discrete breathers can exist in anisotropic spin chains with additional onsite anisotropy [17[ ■ 
In any case, the resultant equations of motion describe an extremely interesting class of dis- 
crete nonlinear dynamical systems and exploration of the underlying dynamical properties 
is of interest both from theoretical and applied physics points of view. 
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The one- dimensional Heisenberg ferromagnetic spin system with nearest neighbor ex- 
change interaction has been shown to possess several completely integrable soliton bearing 
systems in its continuum limit (i) the pure isotropic case (ii) the uniaxial anisotropic 

case j2c|], and (iii) the biaxial anisotropic case These systems also have a strong connec- 
tion with the nonlinear Schrodinger equation |22]. However, till date no exactly integrable 
discrete dynamical Heisenberg spin system has been identified in the literature, although 
a variant of the system, namely the Ishimori spin chain, is known to be completely inte- 
grable 231]. It is generally expected that the discrete anisotropic Heisenberg spin chain is 
a nonintegrable nonlinear dynamical system. Yet, as we show in this article, a number of 
interesting exact periodic and stationary structures, including domain wall type structures, 
for the fully anisotropic system (XYZ case as well as the limitng XYY and planar XY cases) 
can be obtained and their properties analyzed using standard techniques. In fact, Roberts 
and Thompson 2J] and Granovskii and Zhedanov in a series of papers 25|, |26|, |27|, have 
obtained special classes of solutions for the anisotropic spin system. In particular, the latter 
authors have shown that the time-independent case of the XYZ anisotropic spin chain is an 
integrable map by relating it to a Neumann type discrete system 26|], see also Ref. 28 1. 

In this paper, by parameterizing the unit spin vector in terms of the basic Lame poly- 

j, 

nomials of lower order |2^ or their derivatives and by a judicious use of various addition 



30[ | we point out that 



theorems and summation relations obeyed by Jacobi elliptic functions 
several classes of explicit dynamical and static structures can be obtained. In the limiting 
cases we obtain linear spin wave solutions and different nonlinear domain wall type solutions 
in a natural way. We study the physical implications of these solutions like the energy spec- 
trum, effect of discreteness such as the Peierls-Nabarro barrier 3l|, |32|, |33|] , linear stability 
and so on. 

The plan of the paper is as follows. In Sec. II, we introduce the dynamical equations of 
motion and introduce certain natural parametrizations of the unit spin vector. In Sec. Ill, 
we obtain two classes of periodic solutions, investigate the associated dispersion relations 
and energy expressions and indicate a semiclassical quantization of these solutions. In Sec. 
IV we report various classes of static solutions for the XYZ, XYY and XY planar models. In 
Sec. V, we obtain the total energy expressions associated with the various static solutions 
and discuss the effect of discreteness including the Peierls-Nabarro potential barrier. In 
Sec. VI, the isotropic case is considered, while in Sec. VII the linear stability of both time 
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periodic and static solutions is investigated. Then in Sec. VIII, we present some explicit 
time dependent solutions for the case when the onsite anisotropy or an external magnetic 
field is introduced. Finally, in Sec. IX we summarize our results. In the Appendix A 
we include some of the relevant addition theorems and summation relations obeyed by the 
Jacobi elliptic functions required for our analysis, while in Appendix B some details on 
semiclassical quantization are given. 



II. THE HEISENBERG ANISOTROPIC SPIN CHAIN 



A. Equation of Motion 

We consider a one dimensional anisotropic Heisenberg ferromagnetic spin chain with the 
spin components Sn = {S^, S^, S^), satisfying the constraint of unit length 

is:r + {sif + (^:)^ = 1, (1) 

modeled by the Hamiltonian 

H = — y^^i^AS^Sj^^i + BS^S^^i + CS^S^^i) — D ^(5*^)^ — H ■ ^ Sn-, (2) 

{n} n n 

where the sum is over the nearest neighbors. A, B and C are the (exchange) anisotropy 
parameters, D is the onsite anisotropy parameter and Ti = (7i, 0, 0) is the external magnetic 
field along the x direction (for convenience). For the XYZ model, A ^ B ^ C, D = and 
for the XY model C = and D = 0. Using the spin Poisson bracket relation [ij] 

{S^, S^}PB = Sijtap^S], a, f3, = 1, 2, 3, (3) 

where 6ij is the Kronecker delta and ea/^^y is the Levi-Civita tensor, for any two functions A 
and B of spins one has 

{A,B}pB = J2 E^^/^^I^S^^^' (4) 

a,/3,7 i=l » '^'^i 

and the equation of motion becomes 

"^^^ Sn X [A{s:^, + s:_,)U B{si^, + si_,)]+ c{s:^, + s:^,)k 

+2DS'J]+Snxn, n = l,2,...,N, (5) 



dt 



where i, j, k form a triad of Cartesian unit vectors. Explicitly, in component form the above 
equation reads 

= CS'^{S^j^i + — BS^^{S^^i + S'^_i) — 2DS'^S^, (6) 



= AS^{S^_^_i + S^_i) — CS^{S^_^_i + 5'^_i) + 2DS^S^ + HS^, (7) 
dS^ 

- = BS^{S^^^i + S^_i) — AS^{S^j^i + S^_i) — HS^. (8) 



dt 

Equations (5) or (6)- (8) can also be obtained as the limiting case of the corresponding 
quantum dynamical equation of motion for the spin operators when h ^ or S oo. In 
either case, the dynamics is obtained by solving the initial value problem of the system of 
coupled nonlinear ordinary differential equations (5) or (6)-(8) along with the constraint (1) 
on the spin vectors, subject to appropriate boundary conditions like Sn ^n-^oo (±1, 0, 0) or 
Sn —^n-^oo (0,0,±1). However, it appears that the system of differential equations (6)-(8) 
is in general nonintegrable. Even then one can obtain several special classes of solutions 
of physical interest by making use of the properties of (Jacobian) elliptic functions and 
parametrizing the spin vector to satisfy the unit length condition (1). As noted in the 



Introduction, some of these solutions were reported earlier by Roberts and Thompson |24| . 
and by Granovskii and Zhedanov [27], which are to be discussed in the following sections; 
however, as we point out in this paper a much larger class of explicit exact solutions can 
be found in a rather transparent manner through appropriate parametrizations of the spin 
vectors. 

Before dwelling upon the discrete chain, it is also of interest to note as pointed out in 
the Introduction that the long wavelength/low temperature continuum limit of Eqs. (6)-(8), 
when the lattice parameter a ^ 0, takes the form (in the D = limit) 

where H = (7^,0,0), JS = ASj+ BSyj + CS,k, + + Si = 1 and A, B, and C 
are the anisotropy parameters. The isotropic case A = B = C is a completely integrable 



soliton system and is equivalent to a nonlinear Schrodinger equation in a geometrical 18| and 



gauge equivalence sense [l9| . So are the uniaxial anisotropic spin chain {A = B C) in the 
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presence of a longitudinal magnetic field [20| and the biaxial anisotropic spin chain without 



the magnetic field 2l|] integrable soliton systems. In spite of the existence of these integrable 
continuum spin systems, the discrete chain (5) remains as a rather difficult problem to 
analyze. 

B. The Parametrization of the Unit Spin Vector 

One way to proceed with the analysis is to start with an appropriate parametrization of 
the unit sphere of spin given by Eq. (1). Obviously natural parametrizations are in terms 
of elliptic functions. For this purpose one can start with the eigenfunctions of the Lame 
equation 

'L}t^ + [E- n{n + l)khn\u, k)U{u) = (10) 

for positive integer n, which are given in terms of Lame polynomials. The lowest order 
(n = 1) polynomials are [3] 

^/'ii oc sn(M, fc), ■?/'i2 cx cn(u, fc), ipi^ (x (hi{u^k)^ (11) 

while the next order ones are (n = 2) 

ip2i oc sn('U, A;)cn('U, k), ip22 oc cn{u, k)dn{u, k), ip23 oc sn{u, k)dn{u, k), (12) 

and so on. Here sn{u,k), cn{u,k) and dn{u,k) are the standard Jacobian elliptic func- 



tions SOj characterized by the modulus parameter k [see also Appendix A for the relevant 
properties of the Jacobian elliptic functions]. Consequently we can choose, for example, an 
appropriate set of parametrization for the unit spin vectors as 

S^ = asn{u,k), Sy = (3 cn{u,k), 5^ = 7dn(u,fc). (13) 

where and 7 are constant parameters to be fixed. The requirement that condition (1) 
should be satisfied requires 

^2 = 1 _ ^2 ^ ^ ^ _ ^2^/2^ = 1 - f , (14) 

where 7 is a free parameter (0 < 7 < 1) and k is the modulus parameter and k' = 



y/1 — k"^ is the complementary modulus. One can easily check that a parametrization 
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S*^ = asn^u, k)crL{u, k), = Pcn{u, k)dn{u, k) and = 7sn('u, /c)dn('u, fc) does not sat- 
isfy the condition (1) for any set of real values of a, [3 and 7. So one can proceed to higher 



order Lame polynomials [29|| for other possible parametrizations. 

One can even proceed with more general parametrizations in terms of two variables such 

as 

5*^ = cn(M, ki), = sn(M, /ci)cn(f , ^2), 5"^ = sn(M, A;i)sn(t;, ^2), (15) 
with two different moduli ki and k2 or even more general forms such as 
_ acn(n, ki) asn{u, ki)sn{v, k2) 



1 — 7sn(M, /ci)sn(t>, A;2) ' 1 — 7sn(M, /ci)sn(t>, ^2) ' 

nz sn(M, ki)sn{v, k2) - 7 ^ n .^.^ 

= l-7snKfcOsn(.,fc2)' " = V ^ " ^ ' ^l^) 

both of which satisfy condition (1). We will also make use of these parametrizations in our 
analysis. 



III. MOVING SOLUTIONS: ANISOTROPIC CASE (D = 0, W = 0) 

We now look for time dependent moving solutions of Eqs. (6) - (8) when the onsite 
anisotropy and magnetic field are absent (-D = and = 0) in the form (13) and (14) with 
the substitution u = pn — ut + 6, so that 

= a sn(pn — ut + 6, k), = P cn(pn — ut + 6, k), = 7 dn{pn — ut + 6, k), (17) 

along with the relations (14). Here p and uj are the wave vector and angular frequency, 
respectively, which are to be fixed in conjunction with Eqs. (6) - (8) and 5 is a phase 
constant. On substituting the expressions (13) for the components of the spin vector Sn{t) 
into the Eqs. (6) - (8) and making use of the addition theorems for the Jacobian elliptic 
functions (see Appendix A), one requires the following conditions to be satisfied: 

2(3-f[Cdn{p,k)-Bcn{p,k)] 
1 — k'^sn'^{u, k)sio.'^{p, k) ' 

^ 2a7dn(p, k)[Acn{p, k) - C] 
1 — fc2sn2(u, fc)sn^(p, fc) ' 

2 2a/3cn(p,fc)[^- Adn(p,fc)] 

u-ik =— — — — , u = pn-ut + 6. (20) 

1 — /c^sn^(ti, /cjsn^(p, k) 
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In each of the above expressions, we note that the variable u occurs exphcitly on the right 
hand sides. Consequently the u-dependent terms in the above relations can be avoided iff p 
and LJ are chosen in one of the following three ways: 

1. = 0, p ^ 0, 

2. modulus parameter k = (linear spin wave solution, see below), 

3. p = 4:K{k), or p = 2K{k) (or integral multiples of the right hand sides), where K{k) 



is the complete elliptic integral of the first kind 301]. 



In the above, case (1), u = 0, corresponds to the existence of static solutions. These are 
discussed in detail in sec. IV below. 

Case (2), k = 0, corresponds to hnear spin wave (or magnon) solutions. More details are 
given in the subsection IIIc below. 

For case (3), note that sn{2K,k) = = sn{AK,k), cn{2K,k) = -1, dn{2K,k) = 1 = 
cn(4i^. A;) = dn(4i^. A;). Further, since | < K{k) < cx) as < A; < 1, p in Eq. ([21]) is 
bounded below by tt or 27i as the case may be. Correspondingly, we can have two families 
of periodic solutions, each of which we will consider separately: 

p = AK{k), or p = 2K{k), (21) 

A. Spatially homogeneous time-dependent solutions 

For the choice p = 4:K{k), the conditions (18) - (20) reduce to 

uja = -2(3-f{C - B), uj(3 = 2a-f{A-C), oo-fk^ = 2a(3{B - A). (22) 

Solving (22), we obtain 

uj = 2W{B-C){A-C), k'= ^~f Ta'tV iB>A>C), (23) 

7 [A - C) 

where 7 is a free parameter (0 < 7 < 1). The corresponding solutions are 



= v^l - 72A;'2 sn{AKn -iut + 5,k) = -^/l--f^k'^ sn{ujt + 5,k), (24) 

SI = v^l -72 cn{4Kn -ujt + 6,k) = x/T^T^ cn(c<;t + 6,k), (25) 
5^ = 7 dn{AKn -ujt + 6,k)=-f dn{ujt + 6,k). (26) 



The above spatially homogeneous and time periodic solution is nothing but the description 



of Poinsot's motion of a rigid body pointed out by Roberts and Thompson [2J]. Each of the 
spins in the lattice precesses about one of the axes with nutation in the same manner. This 
is depicted schematically in Fig. 1, where all the spins precess parallel to each other. Note 
that the axes stand for the three spin components S^, S"^ and S^. 




FIG. 1: Spatially homogeneous time-dependent solution for the spin vectors. 

The total energy associated with the spin precession of nearest neighbor spins in a 
periodic lattice can be evaluated by substituting the solutions (24) - (26) into the energy 
expression (2) (with D = and = 0) and making use of the summation relations of the 
Jacobi elliptic functions (Appendix A): 

N 

E = — 'y~^JAS„{S„_^_i + S^_i) + BS^{S^^i + 5'^_i) + CS^{Sn^i + S^_i)] 

n=l 

- -2N{BP^ + = -N[B + {C - B)-i'^] 

= -2N[B - {B - C)-i% {B>C, 0<7<1) (27) 
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so that the energy per site becomes 

e=^ = -2[B-{B-Cm- (28) 

Further, since each of these spins evolves identically, the spin chain may be treated to be 
equivalent to independent rigid bodies executing synchronous periodic motions of the form 
(24) - (26). Consequently, each of the above spin motion can be quantized semi-classically 
using the Bohr- Sommerf eld quantization condition 

V^dq^= (^i + ]^h, rii = 0,1,2,..., i = l,2,...,Ar (29) 

where the canonically conjugate variables 

Pn = Sl g„ = arctan ( H ) , n = 1, 2, iV. (30) 



A.. 

Using the explicit forms for S*^, 5*^ and given in (24) - (26), g„'s and p„'s may be expressed 
in terms of elliptic functions (see Appendix B for details). Carrying out the integral over a 
cycle of period AK, one obtains the following transcendental equation for the quantization 
of the amplitude 7: 




n(-^,fc)-(l-f)^(A:) 



n, = 0,l,2,..., ^ = l,2,...,iV. (31) 



Here n(x, k) is the complete elliptic integral of the third kind 30|. Equation (31) is a tran- 
scendental equation in the amplitude parameter 7. For each value of the quantum number 
Hi (=0, 1, 2, ...), the solution 7^. can be found by solving numerically the transcendental 
equation (31). It may be noted that such a semiclassical quantization procedure by solving 
transcendental equations involving all the three complete elliptic integrals has been carried 
out successfully for isotropic anharmonic oscillators with two and three degrees of freedom 
3^ and for the two center Coulomb problem jssl. Then using the resultant allowed set of 
values of the amplitude {7ni}, r;,j=0,l,2,... , i = 1, 2, A^, in the classical energy expression 
per site (28), the corresponding quantized energy spectrum can be evaluated. Consequently, 
the full spectrum associated with the solutions (24) - (26) of the total lattice can be eval- 
uated by associating quantum numbers as {ni, ?t,2, nTv} with the full lattice. Complete 
details will be published elsewhere. 
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B. Spatially oscillatory time periodic solutions 



Now taking the possibility p — 2K{k) in Eq. (18) for the wave vector, and using it in the 
conditions (19) and (20), we obtain the relations connecting the unknowns a;, 7 and k as 

ua^ -2/3'y{B + C), uj(3 = -2a-f{A + C), u'yk'^ = -2a(3{B - A). (32) 

Solving these equations, we obtain 

u;^2W{A + C)iB + C), = ^ (f^) ' P = 2XW- (33) 
The corresponding spatially alternating time periodic solutions are 

= yi - -f^k'^ sn{2Kn -ut + 5,k)^ (-1)"+Vl - 7^^'^ sn(u;t + S,k), (34) 



Sy = cn(2i^n -u;t + S,k)^ (-1)"^! - 7^ cn{ujt + 5,k), (35) 

-S^ = 7 dn{2Kn -ujt + S,k)^^ dn{ujt + S,k). (36) 

The solution (34) - (36) is depicted schematically in Fig. 2. Note that the x and y compo- 
nents of the alternate spins flip and next nearest neighbors evolve in parallel. Also, these 
solutions have no counterpart in the continuum limit of the lattice. 

Again from the energy expression in Eq. (2), using the above solution (34) - (36), we 
obtain the total energy of the system for a periodic lattice of N spins (with N even) 

E^N[B-{B + C)-f% (37) 

Correspondingly the energy per lattice site is 

e^[B-(B + Ch% (38) 

which is greater than the energy of the uniform periodic solution, see Eq. (28). Consequently, 
the present solution constitutes an excited state of the system. 

Since every other spin evolves identically, a lattice of N spins {N even) may be split 
into two sublattices of A^/2 spins, and each member of the sublattice evolves identically. 
Consequently, the semiclassical quantization condition may be given separately for each 
member of the two sublattices: 

j> Pi,idqi,i^ {^i,i + ^h, = 0,1,2, i = l,2,...y 

n2,i + ^h, 712,^ = 0,1,2,..., i = l,2,...y (39) 
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j> P2,idq2,i = ^ 



z 




X 



FIG. 2: Spatially oscillatory time periodic solution. 

and the energy expression may be correspondingly quantized with qi and pi chosen in the 
form (30). From (34) - (36) one can choose for the first sublattice the spin solutions cor- 
responding to n odd and for the second sublattice corresponding to n even in the sohition 
(34) - (36). Then evaluating the integrals (39), which result in expressions essentially of the 
form (31), the energy expression (37) can be quantized, as discussed earlier. 

C. The linear and nonlinear magnon solutions and dispersion relations 

In the uniaxial anisotropic case A — B < C, from the expression (33) for k^, we find that 
k = 0. Consequently, we have the magnon solution. Considering the case (2), /c^ = in 
Eqs (18) - (20), we find that here one has the standard dispersion relation for the uniaxial 
anisotropic case {A = B < C) 

CO = 2-f{C - Acosp), A = B<C, 



12 



corresponding to the linear magnon solution 

= ^1-72 sm{pn -ujt + 6), (40) 
S'^ — a/I — 7^ cos(pn — ut + S), 

which was noted in ref. [18]. 

It is also now instructive to analyze the nature of dispersion relations (33) underlying the 
nonlinear magnon or elliptic function propagating spin wave solutions (34) - (36). 

First we note that these solutions (given by (34) - (36)) in the limit ^ 0, reduce to 
the above linear spin wave solutions (40) with the specific value p = 71 and uj = 2'~f{A + 
C). However, localized solitary wave solutions (for k = 1) do not appear in the moving 
case because p — > 00 (In fact, in the limit k ^ 1, one gets the ground state solution 
(±1,0,0). Prom the expression (33), we note that for 7 = 1, A; = and for 7 = 7^j„ = 
y^{B — A)/{B + C) one gets A; = 1. In other words, 1 > 7 > 7mm- Defining ( — (B — 
A)/(A + C) with B> A implies 



This leads to the dispersion relation [see Eq. (33)] 



{B-A){B + C) 
^ ^ cTP ' ^^^^ 

with p = 2K{k). Note that tt < p < 00 when < A; < 1. For A; ~ 0, K{k) = (7r/2)(l + A;V4) 
and we get the dispersion relation for the magnons (40) as 



(B-A)(B + C) 
^-hl (C-4)+4p/. ■ 

which is finite at p = tt and zero as p — > 00. Similarly, for A; ~ 1 (but not equal to one), we 
have K{k) — ln{A/k') and we get the dispersion for the soliton hke structure as 



^ g / {B-A){B + C) 

"^-'\/a + C)-16exp(-p)' (^^^ 
which is finite at both p — tt and as p — > 00. The above dispersion relations in the limiting 
cases of the parameter k in the allowed region of p, tt < p < 00, are plotted in Fig. 3 for 
a simple choice of the anisotropy parameters. Dispersion curves for other values of k lie 
between these two limiting curves. 
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FIG. 3: Nonlinear spin wave dispersion relations in the limiting cases (i) A; 0: thick line and (ii) 
A; 1: thin line, for the parameter values A = 2.45, B = 2.65 and C = 0. 

IV. STATIC SOLUTIONS 

Considering Eqs. (18) - (20), one can easily note that the term [1 — /c^sn^(-u, k)sn'^{p, k)] 
in the denominator on the right hand sides is "harmless" provided a; = [case(l)] on the 
left hand sides. This actually corresponds to the static case (or time independent case) of 
the equation of motion (6) - (8) for D = and H — 0. In this case the conditions (18) - 
(20) become 



Then there are three possibilities. 

A. Nonplanar static structures for XYZ and XYY models 

In this case S^, S^, ^ 0, which implies that a, /3, 7 7^ 0. The static, periodic solutions 
of such a discrete Heisenberg chain are then obtained from (13) and (14) with u—^m^^ as 



/37[Cdn(p, k) - Bcn{p, k)] 
'ya[Acn{p, k) - C)] 
a(3[B - >ldn(p, k)] 



0. 



(45) 
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(46) 



with the fixed modulus k and fixed wave vector p given by 

yl2_c'2 ' dn(p,A;) = -. (47) 

In (46), 5 is a constant phase factor. Note that the above expressions (47) follow from Eq. 
(45) for a 7^ /5 7^ 0. Here p denotes a wave vector. In the special case of an XYY model, i.e. 
A ^ B = C , which implies k = 1 [from Eq. (47)], we get the localized single soliton (kink- 
and pulse-like) solutions or domain wall structures: 

5^ =tanh(pn + 5), Sl = ^l- 72 sech(pn + 5) 5^ = 7 sech(pn + 5), (48) 

with sech(]9) = B/A. The above domain wall structure is depicted schematically in Fig. (4a). 
Alternatively, for the XXY model, A = B ^ C (i.e. /c = 0), we get the linear excitations 
(frozen magnons): 

^^• = V^T^ sin(pn + 5), = v^T^ cos(pn + 5), ^^ = 7. (49) 



B. Planar (XY) case 



There are two other possibilities from Eq. (45) corresponding to the planar (XY) case. 
Case (i): In this case a = (3 = 1 and 7 = which implies S^ = Q and S^, S-{ 7^ 0. Again 
from Eqs. (13) and (14) the solution is given by 



= sn{pn + 6,k), = cn{pn + 6,k), 



(50) 



provided dn{p, k) = B/A, where the modulus k {0 < k < 1) and the constant 6 are arbitrary, 
which follows from Eq. (45). As in the XYZ case one can obtain domain wall structure in the 



XY case also by taking the limit k —>■ 1, namely 5*^ = tanh(pn + 6) anc 
The above solution was already reported by Roberts and Thompson 



= sech.{pn + 6). 



and Zhedanov 



27| 



24j and Granovskii 
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Case (ii): In this case a = k, P = and 7 = 1 which imphes S-{ = and S!^, 7^ 0. 
The solution is now given by 

S^ = ksn{pn + 5,k), = dn{pn + S, k), (51) 

provided cn(p, A;) = C/A. Here also the modulus < A; < 1 and the constant 5 are 
arbitrary. The domain structure in the k ^ 1 limit is again given by 5*^ = tanh(pn + 5) and 
= sech(pn + 6). This solution was also reported in references 24 1 and 



C. Another class of Nonplanar XYY structures 

Now making use of the more general parametrization (15), one can easily check that in 
the XYY case, that is B = C ^ A with the second variable v fixed as a constant, one can 
identify the following static solution: 



S^ = cn{pn + 6,k), = 'j sn{pn + 6, k), = y^l - -f^ sn{pn + 6, k), (52) 

provided dn{p,k) = A/B, while the modulus parameter k {0 < k < 1), and the constants 
7 = cn(f , ^2) and 6 are arbitrary. In the limit k ^ 1, one obtains the domain wall structure 

= sech(pn + 6), 5^ = 7 tanh(pn + 5), = ^1 - 7^ tanh(pn + 5), (53) 

with tanh(p) = A/B. The above type of domain wall structure is depicted schematically 
in Fig. (4b). It is also of interest to note that none of the above static solutions survive in 
the continuum limit of the lattice and they are all patently structures belonging to discrete 
lattices. 

D. Integrability of the static case 

n 

In their important work ^], Granovskii and Zhedanov have shown that the static case 
of Eq. (5), namely 

X JiSr,+l + Sn-l) = 0, (JSr,) = AS^ + + CS^^, (54) 

or equivalently 

Sn+l + = ^nJ ^Sn, (55) 
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FIG. 4: Localized static single soliton/domain wall structures in the nonplanar case of XYY model 
with 7 = 0.5: (a) Spin structure given by Eq. (48), (b) Spin structure given by Eq. (53). 

where the Lagrange multiplier 

An = ^ > (56) 

is equivalent to a discretized version of the Schrodinger equation with two-level Bargmann 
type potential or a discrete analog of a Neumann system 28|. Some explicit solutions have 
also been reported in their works. 
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V. ENERGY VALUES OF STATIC STRUCTURES AND PEIERLS-NABARRO 
POTENTIAL BARRIER 



For each of the static structures discussed in Sec. IV, the total energy can be obtained 
exphcitly by making use of the various summation formulas given in the Appendix. In the 
following we indicate how the total energy for the lattice can be obtained for the nonplanar 
static structure (46) and then present the final results only for the other cases. 



A. Total energy of the nonplanar static structure (46) 

For the nonplanar static structure (46), the total energy 
E = — ''^[AS^S^^i + BS^S^^i + CS^S^_^_i] 

{n} 

= — - ^^{y4a^sn(u, k) [sn(u + p,k) + sn(u — p, k)] + B0^cn{u, k) [cn{u + p, k) 



+cn('U — p, k)] + C7^dn('U, A;)[dn('U + p, /c) + dn(M — p, k)]}. 



(57) 



Now using the identities derived recently by Khare et al. 
functions, summarized in Appendix A, Eq. (57) can be rewritten as 

N N 



36| for the products of elliptic 



iV iV r 

X [Z{p{n + 1) + 5, A;) - Z{pn + 5,k)- Z{p, k)]. 
Since the elliptic zeta function [SOi] satisfies 

N 

Y,[Z{p{n + l) + 6,k)- Z{pn + 6, k)] = 0, 



BpMn{p, k) C^\n{p, k) 
sn(p, k) 



(5J 



(59) 



n=l 



see Ref. [361], Eq. (58) reduces to 



E = -N\B(3\n{p, A;)+C7Mn(p, k)]-N-^^^^[-Aa^+B(3Mn{p, k)+C^\n{p, k)], (60) 

k^sn[p, kj 

which is independent of the constant phase factor 6. Consequently, the Peierls-Nabarro (P- 
N) potential barrier vanishes. Since a = a/1 — (3 = a/1 — 7^, k"^ = (A^— -B^)/(A^— C^) 

and dn(p, k) = B/A, see Eqs. (46) - (47), the energy expression (60) can be rewritten as 

:b^ - c^] 



E = -N^ - N^^^^Zip, k) 



A 



{A^-B^y 



(A2 - C2) 



(61) 
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Note that in the above expression 7 is a free parameter, while all the other quantitites are 
fixed and the energy is a quadratic function of the free parameter 7. 

In general, this energy may depend on the location of the soliton, i.e. 6. There is an 
energy cost associated with moving a localized mode, e.g. a soliton or breather, by half a 
lattice constant in a discrete lattice. Alternatively, there is a periodic dependence of the 
energy of a soliton on its position with respect to the lattice sites. This is called the Peierls- 



Nabarro barrier 
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32l |. The effects of discreteness such as the P-N barrier and the spin 



barrier (i.e., the total spin, which is an integral of motion, depends on the location of the 
soliton) may be studied from the total energy expressions [33[ . In the present case the P-N 
barrier vanishes as the energy expression is independent of the location of the soliton 6. 



B. Total energy of other static structures 



In the following, we only give the final form of the total energy expressions for the other 
static structures. 

(i) Planar XY case: 

(a) Structure (50): 

E = -NBcn(p, k) + , f ^\ , [A - Bdn(p, k)]. (62) 

k^sn{p,k) 

Here dn(p, k) = B/A, while the modulus parameter k is arbitrary. Consequently the energy 
expression (62) can be expressed as 



where (0 < A; < 1) is the free parameter, 
(b) Structure (51): 

E = -NCdnip, k) + ^^^^[Ae - Ccnip, k)]. (64) 

Here cn{p, k) = C/A, while k is arbitrary. Then the expression (64) can be written as 

E = -N^^A^ - k\A^ - C^) +N ^^ ]^' " (65) 

(ii) Nonplanar XYY structure: For the solution (52), the total energy expression can 
be deduced as 



E = -NAcn{p, k) +N 



2dn{p,k) ^^_2 , r^n .,2m 1 



Ak' -[BY + Cil-Y)] 



sn{p, k) k'^sn{p, k) 



Z{p,k). (66) 
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In Eq. (66), dn(p, k) = A/B while k is arbitrary. Then we have 



E 



-N 



A _ 



B 



k 



- 



(67) 



Note that in the above equation both k and 7 are free parameters. All the above energy 
expressions are independent of the phase constant 6 and so the P-N potential barrier in 
these cases is also absent. 



VI. ISOTROPIC MODEL {A = B = C = I) 

No moving/time dependent nonlinear structures can be found in this case except for the 
spin wave (magnon) solutions (see below). This can be easily checked by looking at the 
conditions that must hold to satisfy Eqs. (18) - (20) which are for k ^ 0, (for k = 0, see 
below) 

ua = -4/?7, ll!(3 = -Aa'j, u'jk'^ = 0. (68) 
Then the only possible structures are the static structures of the following form. 



A. Planar model: static solutions 



In the special case of an isotropic planar model {A = B = C; = 0) the limiting elliptic 
function solutions for 7 = (and u = 0) are 

S^ = sn{2Kn + d,k) , SI = cn{2Kn + 5,k) , ^^ = 0. (69) 

The modulus parameter k and the phase constant 5 are arbitrary. These solutions were 
obtained previously [24] ■ However, the static solutions analogous to Eq. (50) do not exist 
in the isotropic case as K{k) — 00 as A; — 1 and the solution (69) reduces to the uniform 
solution Sn = (±1, 0, 0). However, if instead /3 = 0, the solutions are given by 

Sl = ksYi {2Kn + 6,k), Sl = 0, S'^ = dn {2Kn + 6,k). (70) 



B. Nonplanar model: propagating solutions 

For k = 0, u = 27(1 — cosp) = 47sin^(|), the propagating linear excitations (i.e. 
magnons) are given by 

S: = sm{pn - cut), SI = cos{pn - cjt), S^ = 7. (71) 
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VII. LINEAR STABILITY 



Next, we consider the stability of both the time periodic solutions, discussed in Sec. Ill, 
and the static solutions in terms of Jacobi elliptic functions, obtained in Sec. IV^. Linear 
stabihty of the time periodic solutions in Eqs. (24 — 26) and Eqs. (34 — 36) is studied 
using the period map. For the homogeneous time periodic solution, it suffices to study the 
stability of spin at any one site. For the spatially oscillatory solution in Eqs. (34 — 36), we 
individually study the stability of two spin vectors, one at an odd and another at an even 
site. Rewriting the spin equation (5) using the complex stereographic variable 



we get 



"dT" Hi + |a.+iP n-|a„_,|V "'I 'A i + |a.+iP i + |n„-il' ) ^ ' 

where Cln denotes the complex conjugate of the stereographic variable. After linearizing 
using the expansion 

= ^on + Sn, (74) 
around the time periodic solution Qon{t), we compute the Floquet matrix Xi such that 

sCt„{T) ) \sn„(o) j 

Here, T = 27i/u is the inherent time period in the two sets of solutions, Eqs. (24) - (26) and 
(34) - (36). If 7„ is the eigenvalue of M., then the solution is unstable if |7„| > 1. Figure 
5 shows the instability regions in the {k — A) plane for the homogeneous and spatially 
oscillatory time periodic solutions. Here k is the modulus of the Jacobian elliptic function. 

In order to numerically study the stabihty of the static solutions of Sec. IV, we perturb 
the solution S° in Eqs. (46) - (47) by a small amount 6Sn{t) such that S„ • dSn = 0, 
\SSn\ <C 1. Upon substituting the perturbed vector 

= S° + SSn (76) 



^ The work in this section was carried out in collaboration with S. Murugesh 
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FIG. 5: Instability regions (shaded regions) in the {k — A) plane for the time periodic solutions 
obtained using the period map. Instability region for the (a) homogeneous time periodic solution, 
Eqs. (24) - (26), with B = 0.9 and C = 0.01, B > A > C, and (b) for the spatially osciUatory 
solution, Eqs. (34) - (36), for the same parametric values. The instability diagram is identical for 
spins at both odd and even sites. 



in Eq. (5), the time evolution is computed numerically. As an illustration, it is found that 
the static solution (46) - (47) is indeed stable for long times for small values of the modulus 
parameter k of the Jacobi elliptic function, and that the solution (46) - (47) is less stable 
with increasing value of k, i.e. the instability sets in at earlier times. Figure 6 shows the 
time profile of the static solution (46) - (47) under a small perturbation. Figure 7 depicts 
the initial and final profiles for easy comparison. Fuller details will be presented elsewhere. 

Finally, it is also of interest to investigate whether the time periodic solutions (24) - 
(26) and (34) - (36) are modulationally stable or not. Recently such modulational stability 
analysis has been performed for special solutions of a number of discrete nonlinear dynamical 



systems, including discrete nonlinear Schrodinger equations, see for example refs. |37l. |38|. 
Such an analysis for the time dependent elliptic function solutions (24) - (26) and (34) - 
(36) is being pursued at present and will be reported separately along with the details of 
the linear stability analysis mentioned above. 
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VIII. ANISOTROPIC SPIN CHAIN IN THE PRESENCE OF ON-SITE 
ANISOTROPY AND A CONSTANT EXTERNAL MAGNETIC FIELD 



Finally, in this section we wish to point out that explicit spin solutions can be constructed 
for even more general situations (i) with additional on-site anisotropy and (ii) external 
constant magnetic field. Brief details are as follows. 



A. On-site anisotropy 

Many of the results discussed in the previous sections are valid even in the presence of 
the on-site anisotropy {D ^ 0) in Eq.(5) or Eqs. (6) - (8). In this case, the parametrization 
(17) in terms of Jacobian eUiptic functions leads to the following conditions instead of (18) 
- (20): 

-ua^ 2/37I [^dn(p,fc)-i?cn(p,fc)] _\ 
f dn(p,/c)[Acn(p,fc)-C] ) 

2 2a/3cn(p,A;)[g-^dn(p,A;)] 
^ l-feW(.,fe)sn^(p,fe) ' --Wi-^t^^. (79) 



Unlike the D = case (sec. Ill A), the above relations can be free from the space-time 
variable [u) only for two choices: (i) modulus parameter k — ^ (linear spin wave solution) 
and (ii) p = \.K{k) (spatially homogeneous time-dependent solution) or p = 2K{k) (spatially 
oscillatory time periodic solution). 

(i) For the case k = 0, from Eq. (77c), we have B = A and from Eqs. (77a), (77b) we 
obtain the dispersion relation 

uj = 2'y{C - D - Acos p), A = B < {C - D), (80) 

associated with the spin wave solution (40) for the present case where the on-site 
anisotropy D ^ 0. 

(ii-a) For the choice p — 4:K{k), Eqs. (77) degenerate to 

u;^2W{C-B-D)iA-C + D), ^ J^'"^^ (81) 
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where the form of the associated spatially homogeneous time-dependent solutions are 
given by Eqs. (24) - (26) for the spin components S^, and but with the above 
expressions for cu and k^. 

(ii-b) Similarly for the choice p — 2K{k), Eqs. (77) lead to the expressions 

u;^2WiA + C-D){B + C-D), = (82) 

for the spatially oscillatory time periodic solutions (34) - (36) with the above expres- 
sions for cu and /c^. 

It may be noted from the relation (77) that no static solution (that is a; = 0) is possible 
for general values of k when D ^ 0. 



B. Constant external magnetic field 

Considering now the full anisotropic chain (6) - (8) with the external magnetic field 
also present, we have not succeeded in finding any explicit solution which generalizes the 
nonlinear structures discussed in the previous sections for the anisotropic spin chain in the 
absence of the magnetic field. However, in specific instances classes of exact solutions can be 
obtained. For example, for the XYY spin chain, with B — C ^ A, D — 0, Eq. (5) or Eqs. 
(6) - (8) admit(s) the following exact solutions. In the presence of an external magnetic field 
H — (if J., 0, 0) the relevant contribution to the equation of motion for the anisotropic spin 

— * — * 

chain comes from S'„ x if: 

= Cms:,^, + SU) - S^iSl^^ + Sl_,)l (83) 

dsy 

- ^S:{Sl^, + SU) - CS:{Sl^, + SU) + H^^Sl (84) 

— CS^{S^^i + S^_i) — AS'^{S^^-^ + S^_i) — HxS^. (85) 
These equations have an exact solution with dn(p, k) = C/A and uj = in the form 

5^^' = sn(pn + S,k), 
S"^ = sin(a;t -|- 7)cn(pn + d,k), 

SI = cos{ujt + 7)cn(pn + 5,k). (86) 
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In Eq. (84), the modulus parameter (0 < A; < 1) is arbitrary, while 7 and S are arbitrary 
phase factors. These solutions can be inferred by generalizing the static spin structures of 
the XYY model discussed in Sec. IV. C. Similar solutions can also be written down when 
the magnetic field is along the y or z direction. 
For the present case the total energy becomes 

771 _ A \ ^ QX QX ^ \ ^ { cy cv _L ^ ij \ 

n n n 

= —A sn(pn + S) sn{p{n + 1) + S) 

n 

—C cn(pn + 5) cn(p(n + \) + 5) — H sn(pn + 5) 

n n 

= -NC cn(p. A;) - ,^^^f ' M + C dn(p. A;)] - if,. V sn(pn + 5, k) 
k sn(p, K) 



_^C VC^-A^+ ^ + iV ^^'~^' z(p, k) + H^Yl ^^(^^ + ^' (87) 



A 



k 



The sum in the last term above, namely Yin sn(pn + 5), is represented as era (5) in Ref. 36 1 
by Khare and Sukhatme and is dependent on the location of the soliton 6. This ensures that 
the Peierls-Nabarro potential barrier is present in the anisotropic spin chain in the presence 
of an external magnetic field. 



IX. SUMMARY AND CONCLUSIONS 



Using the summation identities 36|] for Jacobi elliptic functions 30| we obtained sev- 
eral classes of static and propagating exact solutions for the classical, anisotropic Heisen- 
berg chain. In the special case of the isotropic planar model we recovered the previously 



known solutions 2^- We explicitly obtained the nontrivial dispersion relations {u vs. p) for 
the propagating solutions and predicted the contrasting features of magnons and solitons. 
Specifically, as p — 00 the magnon frequency goes to zero whereas the soliton frequency 
reaches a nonzero value. These dispersion relations can be measured via neutron scattering 
in the quasi-one dimensional materials realized by anisotropic Heisenberg chains. It would 
be instructive to explore whether similar exact solutions can be obtained for the analo- 
gous quantum Heisenberg models. The effects of discreteness, e.g. Peierls-Nabarro barrier 
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32, 133(1 and spin barrier 33|] , may be important in anisotropic spin chains. The discrete 



equation of motion is non-integrable in general. However, the static version in the absence 



25 



of an external field is an integrable system 25|. It is instructive to explore semiclassical 



quantization of the 2K versus AK periodic solutions in terms of anharmonic oscillators 



39l |. The solutions expressed in terms of sn(x, fc), cn(x, fc) and dn(x, A;) correspond to the 
n = 1 Lame functions. We have explicitly checked that n = 2 Lame functions do not give 
exact solutions. However, it is conceivable that n = 3 Lame functions [29] may lead to a new 
class of exact solutions. In addition, there may be another class of exact solutions with a de- 
nominator also containing elliptic functions. We are presently exploring these solutions also. 
To conclude, we wish to state that the classical anisotropic Heisenberg spin chain admits 
very interesting static and dynamic structures and more work is needed in this direction to 
identify all of them. 
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APPENDIX A: JACOBIAN ELLIPTIC FUNCTIONS 

Basic elliptic function properties: 

sn^{u,k) + cn^{u,k) = 1, dn^(M, fc) + Psn^(M, A;) = 1. (Al) 

Addition theorems: 

. , 2sn(-u, fc)cn(f , fc)dn(f , fc) , ^ ^, 

sn M + +sn M-t;,A; = \; m 2/ 

1 — A;"'sn"'('U, A;)sn^(f , k) 

/ , N / , N 2cn(u, A;)cn(t>, /c) , . , 

cn(« + ^, fc) + cn(7/ -v,k)= \ \ ^ ;/ (A3) 

1 — K^sn^^-u, A;)sn"^(f , k) 

1 / T\ 1 / 7x 2dn(u,k)dn(v,k) 

dn{u + v,k) + dn{u -v,k) = , / ' ' ' \/ A4 

1 — A;^sn^(M, A;)sn^(u, A;) 

Product relations: 

msn(x, A;)sn(x + a, fc) = — ns(a, k)[Z{x + a) — Z{x) — Z{a)], (A5) 
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mcn(x, /c)cn(x + a,k) = mcn(a, k) + ds(a, k)[Z{x + a) — Z{x) — Z{a)], (A6) 

dn(a;, fc)dn(x + a,k) = dn(a, /c) + cs(a, k)[Z{x + a) — 2'(a;) — ^(a)], (A7) 

where Z{x) = Z{x,k) is the Jacobi or elhptic zeta function, and ns{x,k) = l/sn{x,k), 
ds(x, k) = dn{x, k)/sn{x, k), cs(x, k) = cn{x, k)/sn{x, k). 
Summation relation: 

N 

{Z[l3e{n + l) + 5,k]- Z[n(3e + 5, A;]} = 0. (A8) 

n=l 

Integration formula: 

where K{k) and n(Q;^, k) are the complete elhptic integrals of the first and third kind, 
respectively. 



APPENDIX B: SEMICLASSICAL QUANTIZATION 



For the spatially homogeneous time- dependent solution (24) - (26), the canonically con- 
jugate variables qi and Pi given by Eq. (30) can be expressed in terms of elliptic functions 

as 

Pi = S- = 7dnti ; u = ut + 6, (Bl) 



Qi = arctan ( I = arctan 



1 — 7^ cnu 
1 — 7^A;'^ snu 



(B2) 



Then the left hand side of the semiclassical quantization condition (29) becomes 30|] 

dn\ du 



7^dn^M 



4 /l-72A;'2 
7V 1-7' 



n 



A; - {l--f^)K{k) 



(B3) 



(B4) 



,(1-7^) 

where K{k) and n(z/, /c) are the complete elliptic integrals of the first and third kind, re- 
spectively. 
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time 



FIG. 6: Time evolution of the Sn component of the static solution, Eqs. (46) - (47), under 
perturbation, (a) k = 0.3, (b) k = 0.9 and (c) ^g= 1. As can be noticed, for small values of k, the 
solution tends to be more stable for long periods of time. 
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FIG. 7: Initial (line-points) and final (line) profiles of the static solution (46) - (47) under pertur- 
bation (Eq. (76)) in (a) Figure 6(b) with k = 0.9 and (b) Figure 6(c) with A; = 1. Instabilities 
start to appear at earlier time as k approaches 1. 
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